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Abstract. An external magnetic field induces large relative changes in the entropy 
of one-dimensional quantum spin systems at finite temperatures. This leads to a 
magnetocaloric effect, i.e. a change in temperature during an adiabatic (de)magne- 
tization process. Several examples of one-dimensional spin-1/2 models are studied 
by employing the Jordan- Wigncr transformation and exact diagonalization. During 
an adiabatic (de)magnetization process temperature drops in the vicinity of a field- 
induced zero-temperature quantum phase transition. Comparing different levels of 
frustration, we find that more frustrated systems cool down to lower temperatures. For 
geometrically frustrated spin models a finite entropy survives down to zero temperature 
at certain magnetic fields. This property suggests frustrated quantum spin systems as 
promising alternative refrigerant materials for low-temperature magnetic refrigeration. 



PACS numbers: 75.30.Sg, 75.10.Jm, 75.50.Ee 
1. Introduction 

The magnetocaloric effect, or lieating/cooling of a magnetic system under an adiabatic 
(de) magnetization process, is related to the isothermal variation of the entropy in a 
magnetic field via a simple thermodynamic relation 

dT_\ {dS/dH)T 
OH J, C ' 

where C is the specific heat at a constant magnetic field. A field variation of entropy can 
be found via a temperature variation of the magnetization: [dS/dH)T = {dM/dT)H- 
Standard examples of magnetic materials with a significant magnetocaloric effect include 
paramagnetic salts, which consist of noninteracting spins, and ferromagnets near the 
Curie point, where {dM / dT)^ is, generally, large [1112] • Besides a fundamental interest, 
the magnetocaloric effect has a great importance for the technique of magnetic cooling. 
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which has been employed over the years to reach temperatures in a sub-Kelvin range (3] . 
A certain progress has also been achieved to utilize this technique for room temperature 
refrigeration jUEj- 

Recently, it has been predicted that an enhanced magnetocaloric effect exists in the 
vicinity of finite-field transitions in a class of geometrically frustrated antiferromagnets 
[H]. Local spins in such magnets reside on special types of crystal lattices, with 
pyrochlore, kagome, and garnet lattices being a few well-known examples. Competing 
antiferromagnetic interactions on frustrated lattices lead to a macroscopic degeneracy 
(entropy) of classical ground states. Large variations of the total entropy occur near 
the saturation field He, where the transition to a fully polarized nondegenerate state 
takes place. Above He in the saturated phase, highly frustrated magnets have a flat 
dispersionless low-energy branch of magnons with energy {H — He). Dispersionless 
excitations can be represented in real space as localized quasiparticles. Such localized 
magnon excitations have recently also been studied in frustrated quantum spin models 
[ZIIH] where it was shown that the corresponding many-particle states are not only exact 
eigenstates of the Hamiltonian, but indeed the lowest ones. Hence, the macroscopic zero- 
T entropy of geometrically frustrated magnets at H = He is robust against quantum 
fluctuations (see also [IHl), whereas below He the quantum (or thermal) order by 
disorder effect should lift the classical degeneracy [ZllII]. It is the condensation of a 
macroscopic number of localized magnons near the saturation field, which leads to the 
enhanced magnetocaloric effect in geometrically frustrated magnets jH]. 

The aim of the present work is two-fold. First, we argue that a large magnetocaloric 
effect appears near continuous phase transitions in a magnetic field for general one- 
dimensional (ID) quantum spin systems. This phenomenon has a simple physical origin. 
The excitation spectrum at the quantum critical point is different and, generally, softer 
if compared to spectra above and below the transition. The density of low-energy modes 
is strongly enhanced in spin chains due to the ID van Hove singularity. Therefore, the 
soft modes increase the total magnetic entropy around the transition point and produce 
a sizeable magnetocaloric effect. We present analytic calculations for two exactly solved 
models: an XY antiferromagnetic chain with a magnetic field applied perpendicular 
to the plane and an Ising model in a transverse field. Then, we study by means 
of exact diagonalization of finite lattices a spin-1/2 Heisenberg chain. Secondly, we 
investigate how the magnetocaloric effect in ID systems changes with increasing degree 
of frustration. For this purpose we consider a frustrated J1-J2 spin-1/2 chain and the 
so-called sawtooth spin chain as a ID toy model for the class of geometrically frustrated 
lattices containing the kagome and pyrochlore lattice. 
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2. XY chain in transverse field 

We begin with a simple spin-1/2 XY antiferromagnet in a transverse magnetic field, 
which is described by the Hamiltonian: 

N 
i=l 

With the help of the Jordan- Wigner transformation jT2I this is mapped to a Hamiltonian 
of free spinless fermions: 

Tt; = \-^ekclck , ek = H + Jcosk. (3) 

k 

In high magnetic fields the ground state is the fully polarized spin state, which describes 
the fermion vacuum. As the magnetic field is lowered, the excitation gap for fermions 
decreases and vanishes at the saturation field = J. At the transition point the 
low-energy fermions have a parabolic dispersion e^+g ^ ^Jq^ for |g| -C 1. 
The heat capacity of a gas of free fermions is 

C = ^Y.^kM'^-nk) , (4) 

k 

where = [e^''/-^ + are the occupation numbers. At the saturation field the specific 
heat has a square-root low-temperature dependence 

3 



C/N ^ 



with C(|) ~ 2.612375. 

The magnetization of an XY chain is given by M = — 'Y^^n^, which reduces 
at zero temperature to 

M/N = --- aTccos{H/Hc) (6) 

2 TT 

(compare with eq. (3.5) of JSI)- The magnetization curve M{H) has a standard square- 
root singularity near the saturation field P^1IT3] . The magnetocaloric effect is calculated 

as 



dS 



OH T2 

' ^ k 

At H = He and low temperatures one finds 
/dS\ ciN 



l^'^ekUkil-Uk) . (7) 



Cl 
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Figure 1. Adiabatic demagnetization curves of the XY antiferromagnetic S = 1/2 
chain. 



The low-temperature expressions for the entropy above, at and below the saturation 
field He are given by 

Ah 



S/N 

S/N 
S/N 



-Ah/T 



V2ttT 

^C(i)(l-2-/^) 

TT T 

3 ■ 



T 
2^ 



H>H, 



H=Hr 



H<H, 



(9) 



^2\Ah\J 

where Ah = H — He and \Ah\ ^ T is assumed. At T = in an arbitrary magnetic 
field the system has a unique quantum ground state and, correspondingly, zero total 
entropy. At T > 0, thermal fiuctuations yield different entropy contributions in the 
above three regimes. Variation of the applied magnetic field leads, therefore, to large 
relative changes in the total entropy and to a large magnetocaloric effect. 

The adiabatic demagnetization curves of an XY antiferromagnet can be found 
either from a direct numerical solution of S{H, T) = const or by numerical integration of 
the differential relation (P). The results are presented in Fig.^for Hi = 2Hc and several 
starting temperatures Tj. At H = He the cooling rate {dT/dH)s has the maximal value 
determined by eqs. ((T)), (jSJ and (jH)). The lowest temperature of an adiabatic process 
is reached at H* < He- At low temperatures the difference between these two fields 
becomes, however, very small. Using eqs. Q we find the following analytic expressions 
for a field-dependence of temperature. At {H — H^) ^ T, temperature changes nearly 
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linearly with {H — He) due to a dominant exponential contribution. In this region 

(10) 

The lowest temperature reachable during an adiabatic demagnetization process from Tj 
and Hi > H^ to Hf = He is exponentially small: 

2{Hi — He 



(H, - HeY 

Tf exp 
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T- 

Below the saturation field, temperature decreases as (H^ — H)^^"^ with H H^ 
[{He — H) ^ T]. In this region one finds 

.1^. (12) 

The lowest temperature reachable during adiabatic magnetization from H = Hi < H^ 
to Hf = He is 

Tf ~ . (13) 

3. Heisenberg chain in magnetic field 

We now discuss the effect of S^-S^ interaction on the magnetothermodynamics of the 
phase transition in the spin-1/2 XX Z chain: 

^ = E (^^'^^^1 + + ^StSU,) - HS!] . (14) 

i 

The case A = 1 corresponds to the Heisenberg chain. The interaction of the z- 
components shifts the fermion energies and the saturation field: 

ek = H - JA + Jcosk , He = J{l + A). 

It also induces an interaction between fermions: 

V = JA J^4ci4+iCi+i 

i 

= ^ 5Z J^cosq- cos{k ~k' + q)]cl^^cl,_gCk'Ck . (15) 

k,k',q 

The interaction between spinless fermions vanishes as k,k',q ^ 0. In the vicinity of the 
saturation field only wave-vectors near k ^ tt are important and, therefore, interactions 
play a minor role [T^15pi6j . The main correction to the free energy at low temperatures 
is the Hartree-Fock contribution, which is estimated as 

SFhf = {V) = JNA{n^ - m^) (16) 

with n = (clci) and m = (c|q+i). At the saturation field one finds 



El — cos k T 
es,/T 1 V 27rJ 

k 



El + cos/c 1 > /-JN / 1 /9\ 



k 
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The temperature expansion of the free energy at H 

, r 

m = -c (i) (1 - 2-^/^) 



He is 



Aci 



T2 



:i7) 



_(27rJ)i/2 27rJ 

where the first T^/^ term is a Fermi-gas contribution and the second term is the 
Hartree-Fock correction. The next order T^/^ correction to eq. (fT7|) is determined 
by deviations from a parabohc dispersion in Sk- The interaction gives a neghgible 
contribution in the thermodynamics of a Heisenberg chain at temperatures T < O.IJ. 
The expression (|T7jl agrees with the exact thermodynamic Bethe ansatz solution for a 
Heisenberg chain f21E|, which, however, gave the second term without numerical 
prefactor. 



4. Frustrated Heisenberg chain near saturation 

The density of low-energy states at the saturation field can be increased by frustration. 
As a first example we consider the J1-J2 chain with exchange interaction between the 
first and second neighbours at the point J2 = ^Ji, where the dispersion degenerates 
from the generic quadratic behaviour to a fourth power, giving rise to a fourth-root 
instead of the usual square-root singularity in the T = magnetization curve at the 
saturation field [UlIiniEn] • 

The J1-J2 Heisenberg chain is described by the Hamiltonian 

n = J2 [^iS. ■ S.+i + ■ S,+2 - HSn . (18) 

i 

The fermion energies in the saturated phase are given by 

Ek = H + Ji{cosk - 1) + J2(cos2A; - 1) . (19) 

Frustrated interaction with the second neighbours does not change the saturation field 
value for J2 ^ \Ji- He = 2Ji. The long- wave length expansion of Sk near the minimum 
point yields 

..,,./^-i/.+ (A_2..),^^(f -A),.. 

At the point J2 = ^Ji, the first g^-term vanishes and the dispersion of fermions is 
described by a quartic parabola: e.j,+q ~ H — He + \J'i^ {Ji = J)- This leads to a 
higher density of states above the gap p{uj) ~ uj~^^'^ than for the standard ID van Hove 
singularity p(tu) ~ uj~^^'^. For a free fermion gas with a general dispersion law Sk = Ak^, 
the entropy behaves as 

+ + ^) (1 - ^ • (20) 

For z = A this gives S ~ T^^^. However, for z > 3 the interaction between fermions 
becomes again a relevant perturbation aX. H = He [IH] such that eq. is no 

longer expected to give the exact low-temperature asymptotics. Thermal fluctuations 
tend to reduce frustration by generating a temperature-dependent quadratic term 
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Ek ~ B{T)k'^ + Ak^, which has to be determined self-consistently. In any case, we 
expect that at low temperatures the entropy of the frustrated chain is higher than the 
entropy of a Heisenberg chain. Since above He the entropy has a similar exponential 
behaviour given by eq. the frustrated chain has larger relative variations of the 
entropy and, consequently, should exhibit a larger magnetocaloric effect. 

5. Numerical results at the saturation field 

In this section we compare the analytic results for the S = 1/2 Heisenberg chain to 
numerical results obtained by exact diagonalization of finite rings and check the scaling 
of the entropy for the J1-J2 Heisenberg model. A spectral representation is used for 
the evaluation of the thermodynamic quantities. In particular, the entropy is computed 
using 

n 

where Z = q-^^/t ^.j^g partition function. The eigenvalues En are obtained by full 
diagonalization of the Hamiltonian. 

Since the CPU time for a full diagonalization grows as the third power of the 
dimension, it is important to use symmetries to reduce the dimension. We have used 
conservation of S^, the z-component of the total spin, spin-inversion symmetry for 
= and S't/(2)-symmetry to reconstruct the spectra in the sector = 1. Then 
the biggest dimension arises for 5*^ = 2. The actual size depends on the spatial 
symmetries, namely translations and reflections, which have also been exploited. Full 
diagonalizations have been performed with up to = 20 sites for all three spin-1/2 
systems to be discussed below. The sawtooth chain has lowest spatial symmetry, hence 
the largest dimension occurs here, namely 12 618 for = 2 and N = 20. Since the 
magnetic field H couples to the total magnetization S^, which is a conserved quantity, 
one computation yields the dependence of the entropy S both on temperature T and 
the magnetic field H. 

Fig. 121 presents a log-log plot of the temperature-dependence of the entropy of the 
Heisenberg chain and the J1-J2 chain with Ji = J, J2 = at the saturation field 
He = 2J for = 12, 16, and 20. At low temperatures, the ground state at He is two- 
fold degenerate, i.e. 5* = ln2 for any A^ as T ^ 0. Accordingly, one observes finite-size 
effects in S/N at low temperatures. At high temperatures, all curves recover the full 
entropy S/N = In 2 pa 0.6931 of a spin-1/2 system. Finite-size effects are still negligible 
at intermediate temperatures and the numerical curves are straight lines in the log-log 
representation of Fig. |21 implying a power law for S/N as a function of temperature. 
Next, we discuss these power laws in more detail and argue that the exponents are 
indeed the expected ones. 

The bold solid line in Fig. |21 obtained from eq. (fTTj) is (see also eqs. and (jSDl)) 

S/N = 0.457874 (T/J)^/^ - 0.1473254 (T/J) . (21) 
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Figure 2. Numerical results for the entropy of the spin-1/2 Heisenberg chain (dotted 
lines) and the J1-J2 chain with Ji = J2 = jj (dashed hues) at the saturation field 
He — 2 J. The behaviour for a free fermion gas (FF) is shown by the bold solid line, 
eq. H21|l . and the thin solid line, eq. H2()() with z = 4 and A = ^J. 

This fits the numerical data nicely in a wide temperature range 0.007 ^ T/J < 0.2. 
Deviations from ()21|) for lower temperatures are determined by finite-size effects in 
the numerical data, whereas at higher temperatures further corrections to (j2ip become 
important. Although indeed the agreement at higher temperatures can be improved by 
including the next correction 0(T^/^), we will not pursue this further here. 

The thin sold line in Fig. El shows eq. ()20|) with 2; = 4 and A = appropriate for 
the frustrated chain. Comparison with the numerical data for intermediate temperatures 
confirms the expected power-law behaviour 5* ~ T^/^. However, the prefactor derived 
from eq. (j^n|) for z = 4 is by about 20% too big in contrast to the Heisenberg chain 
{z = 2). This behaviour of the prefactor reflects the fact that interactions are relevant 
for z = 4 and partially reduce the entropy compared to the free-gas expression eq. (j^n|) . 

6. Anisotropic XY chain in transverse field 

Low crystal symmetry of real (quasi) one-dimensional magnetic materials is compatible 
with various types of spin anisotropy. If the magnetic field is applied in a general 
direction, the anisotropy leads to a nontrivial behaviour by breaking the conservation 
of the total magnetization. It is, therefore, important to investigate the effect of spin 
anisotropy on the magnetocaloric effect. The simplest kind of anisotropy is anisotropic 
exchange interaction. We consider the following Hamiltonian 

n = J2 {J^S:St,, + JyS^Sl, - HSl) , (22) 
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Figure 3. Dependence of the magnetocaloric effect on anisotropy. The number near 
each curve indicates the corresponding value of the anisotropy parameter 7. 



which at Jx = Jy reduces to the XY chain and at Jj^ = describes an Ising chain in a 
transverse field. Introducing new parameters J = ^{Jx + Jy) and 7 = {Jx — Jy) / {,Jx + Jy) 
and applying the Jordan- Wigner transformation, the Hamiltonian (j^^ is reduced to 



HN 



i"yJ . 



H = -—^ + y ]{H + J cos k)c\ck + ^^sin/i;(c|,c^^- C-kCk) ■ (23) 

k 

The Bogoliubov transformation from original fermion operators to new operators a^: 

Cfc = MfcOfc + ivk0^_^ , c_fc = Uka^k - ivka\ 

diagonalizes the quadratic Hamiltonian (j23|) under the two conditions: (i) u 
and (ii) 

7 J sin k 

2ukVk -- 



1 



v/(i/ + JcosA;)2 + 72j2 sin k 
The transformed Hamiltonian has a simple diagonal form: 

1 



ek = \{H + Jcosk)^ + -f^J^sm'k. (24) 



The anisotropic XY chain (j22j) has a phase transition at He = J. Near the 
transition point the long-wavelength spectrum is given by {Ah = H — Hf) 



^'Af, + 7VV. 

The dispersion acquires a relativistic form aX H = He- ^Tr+g ~ 'jJq- The expressions for 
the free energy, entropy and the specific heat are the same as for a gas of free fermions. 
The magnetization of the anisotropic XY chain is 

H + Jcosk /I 

~k u 



M 



rik 



(25) 
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At zero temperature Uk = and the total magnetization does not saturate above He- In 
the case of the Ising chain 7 = 1, the magnetization has a universal value M{Hc) = l/vr 
at T = 0. The square- root singularity near He is replaced with a weaker logarithmic 
behaviour: dM/dH ~ \n.{'K'y / /S.h)- 

The expression for the magnetocaloric effect changes compared to the XY chain 

(H) ^ = + ^ - ^'^) • (26) 

Due to a linear dispersion law the density of low-energy states in the anisotropic chain 
remains finite instead of the square-root divergence typical for isotropic spin chains. As 
a result, the entropy at H = He depends linearly on temperature: 

This leads to a reduction of the magnetocaloric effect compared to the isotropic case. 
The cooling rate at the saturation field also changes its temperature dependence for the 
anisotropic chain: 



OH J g 272 
dT \ 2ci 



anisotropic XY , (28) 
isotropic XY . (29) 



OHJs 3C(|)(1-^) 

The dependence of the adiabatic demagnetization process on the anisotropy is shown 
in Fig. ini The lowest temperature reachable at H ^ He grows significantly from the 
isotropic case 7 = to the strongest anisotropy 7 = 1. Nevertheless, even an Ising chain 
exhibits sizeable variations of temperature. Another difference with the isotropic chain 
is a stronger temperature increase below H^. This is related to a gap opening in the 
anisotropic chain, which leads to a significant reduction of entropy. 



7. Universal properties 

Universal low-temperature scaling at second-order quantum phase transitions is related 
to the dispersion of gapless excitations at the transition point. In particular, the 
scaling laws (j28|) . (j29|l remain valid in the general case: {dT/dH)g~^'' = 0(1) for 
t/(l)-symmetry about the field direction and {dT / dH)^^^" = 0{T) for Z2-symmetry 
about the field. In fact, this is not only true close to the saturation field, but also in 
the vicinity of plateaux in the magnetization curve [011211122] ■ Since these field-induced 
transitions are universal (at least in one dimension), the scaling laws given by eqs. ()28p 
and ()29|l should apply to generic second-order transitions in one dimension. 

Our results can be related to a general scaling analysis of the magnetocaloric 
effect close to quantum critical points |2S1- First we note that e.g. the power laws 
in eqs. (1201) and ((211) can be written as C, 5 oc T'^/^ at H = H^ and C, 
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S oc T^^ \H - H^l^^'^'^oz) jj ^ jj^^ ^gj^g notations of The case of U{1)- 
symmetry is described by the exponents {d = 1) 

z = 2, iy=^, (30) 

7/0 = 1 in a gapless phase and ?/o = oo in a massive phase. For the case of Z2-symmetry 
the exponents are given by {d = 1) 

z = 1, 1^ = 1, yo = oo. (31) 

With these identifications, our resuhs (fTUI) and ()12p for the isotropic XY chain are 
consistent with the predictions from general scahng arguments |^|. Furthermore, eq. 
^ obeys the general scaling form {dT/dH)^^"' oc T^-^i^^). 

The situation with only Z2-symmetry is an exception in so far as the general scaling 
arguments [22] do not directly lead to eq. (j2HI)- In this case, the leading term cancels at 
H = He and eq. is obtained from the next order in T. For Z2-symmetry, H ^ He 
and T <^ \H — Hc\, the general scaling predictions [231 amount to 
1 fdT\ _ 1 

T\dH) ,~ir^,' ^^^^ 

which coincides with the result for the f/(l)-symmetric system in the gapped state eq. 
()10|) and agrees with an explicit computation for the anisotropic XY chain. 



8. Sawtooth chain 



The class of models with localized magnon excitations at the saturation field includes 
the spin-1/2 Heisenberg models on kagome and pyrochlore lattices jZllHlini- In this section 
we study the magnetothermodynamics of a one- dimensional representative of this class, 
namely the sawtooth chain (called also a A chain). The sawtooth chain is a chain of 
corner-sharing triangles with two different exchanges for the base-base and base-vertex 
edges, see Fig. HJ For the special ratio of the coupling constants J2 = \Ji the sawtooth 
chain supports localized magnon excitations, giving rise to a plateau in the T = 
magnetization curve at half the saturation field followed by a jump to saturation in the 
spin-1/2 model 0^. Below we always assume the choice of parameters J2 = |<^i with 
Jl = J. 

In the fully saturated phase at high magnetic fields the sawtooth chain has two 
branches of single spin-flip excitations with energies 

eik = H -2J , e2k = H -^{l-cosk) . (33) 

Excitations from the low-energy dispersionless branch correspond in real space to 
localized spin-flips residing in valleys between two triangles, see Fig. E] Such a 
localization is a direct consequence of the frustration imposed by triangular topology 
and resembles to a large extent similar behaviour of magnons in kagome, pyrochlore 
and other frustrated lattices. Localized magnons in the sawtooth chain condense below 
the saturation field H^ = 2J and form at zero temperature a dense magnon crystal 
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Figure 4. The sawtooth chain. LocaUzed magnons, which exist for Ji — 2J2, are 
trapped within one valley between two adjacent triangles with amplitudes indicated 
near each site. 



with one spin-flip occupying every second valley between triangles [HIH]. The magnon 
crystal ground state is two-fold degenerate and breaks the translational symmetry. In 
the following we discuss the low-temperature thermodynamics of the sawtooth chain in 
the vicinity of the saturation field, which is dominated by excitations from the lower 
branch. 

We consider a chain which consists of L triangles (A^ = 2L spins). The lowest energy 
of a state with n excited magnons, n < L/2, is equal to ne , e = H — He = eik- The 
corresponding states are formed by nonoverlapping localized spin-flips, which occupy 
spatially separated valleys. This suggests a hard-dimer representation of the lowest 
energy n-magnon states so that the presence of a dimer on a given lattice site of the 
basal chain in Fig. |3] forbids occupation of two neighbouring sites by other dimers. Such 
a simple construction of the lowest energy states in every ra-magnon sector is complicated 
by the fact that the corresponding states are not orthogonal due to an overlap of the 
wave-functions of spin-flips trapped in two adjacent valleys. The states are, however, 
linearly independent and can be properly orthogonalized. Besides such lowest-energy 
localized states there are also scattering states in every n-magnon sector. The scattering 
states of n magnons are separated from ne by a finite gap A^. For the two-magnon sector 
of the sawtooth chain the variational calculations estimate A2 ~ 0.44 J and suggest 
An = 0(A2) in other sectors j21]. The present numerical calculations confirm this 
picture yielding A2 — 0.42 J, whereas for the magnon crystal state on the 1/2-plateau 
Al/2 — 0.25 J. Neglecting at low temperatures exponentially small corrections from the 
scattering states, we are left with the calculation of the free energy of hard dimers. This 
is a classical problem without dynamics and can be calculated exactly (see pHl25j ). The 
classical partition function of the hard dimers is 



where the index i runs over the sites of the basal chain and cTj = 0, 1 denotes the absence 
or presence of a dimer on a given site. The exclusion principle for nearest-neighbour sites 
is imposed by the last term. The partition function ()34j) can be rewritten as Z = trT^, 




(34) 



Magnetocaloric effect in one- dimensional antiferromagnets 



13 



T being the individual transfer matrix given by 




(35) 



This immediately yields the free energy and the entropy normalized per triangle 



which amounts to 34.7% of the total entropy 2 In 2 of the spin-1/2 sawtooth chain. 
The same zero-temperature entropy of the sawtooth chain at H = He has been 
independently obtained in [T^. It is also known to express the entropy of an Ising 
antiferromagnetic chain at the saturation field [SHI- The entropy Sh=Hc has no 
temperature dependence in the hard-dimer approximation. A temperature dependence 
does arise from neglected multimagnon states and its magnitude can test the validity of 
the hard-dimer representation for the low temperature thermodynamics of the sawtooth 
chain. 

The results on the temperature dependence of the entropy obtained by exact 
diagonalization for = 20 sites and by the hard-dimer mapping are compared in Fig. El 
for several magnetic fields. For H > He = 2, the numerical results are in excellent 
agreement with the hard-dimer approximation at low temperatures consistent with a 
gap A ~ 0.2 J to states neglected in the hard-dimer picture. The magnon crystal arising 
for H < He is two-fold degenerate. Hence, one obtains a non-extensive entropy 5* = ln2 
at T = 0, giving rise to finite-size effects in S/N for H < H^ (see numerical curve 
for = 1.9 in Fig. E}. In this regime, the hard-dimer approach is superior at low 
temperatures since it does not suffer from finite-size effects. The inset of Fig. El shows 
the field dependence of the specific heat at T = 0.05J. For H > H^ the hard-dimer 
approximation agrees very well with the numerical results for = 20. For H < H^, 
both approaches become approximate since the numerical result (full line) is affected 
by finite-size effects and the hard-dimer approximation neglects states that are relevant 
in this regime. The two peaks structure in the specific heat on both sides of the critical 
field obtained in the two approaches is typical for field-transitions with finite entropy at 
H = He^. 

9. Comparison of different degrees of frustration 

In the preceding sections we have considered the spin-1/2 Heisenberg model, the J1-J2 
chain at J2 = | Ji, and the sawtooth chain at J2 = | Ji. These models can be considered 




(36) 



At H = He the total entropy is 




(37) 



Magnetocaloric effect in one- dimensional antiferromagnets 14 




T/J 

Figure 5. The temperature dependence of entropy for the sawtooth chain for different 
magnetic fields. Full lines are numerical results for = 20 sites, dashed lines are from 
the hard-dimer approximation. The number near each line indicates the corresponding 
value of the external field. The inset shows the field dependence of the specific heat C 
a.t T — 0.05 J from numerical (full line) and the hard-dimer (dashed line) calculations. 

as examples for increasing degree of frustration, as is reflected e.g. by an increase of the 
entropy at low temperatures at the saturation field H^. In this section we will compare 
their entropy and adiabatic (de)magnetization curves in more detail. 

Fig. El shows the entropy as a function of temperature and magnetic field for the 
Heisenberg chain, the J1-J2 chain with Ji = J, J2 = \ J ■, and the sawtooth chain with 
Ji = J, J2 = |J. The saturation field is i^c = 2J in all three cases. The contour 
lines show curves of constant entropy, i.e. adiabatic (de) magnetization curves. The 
global behaviour of all three models exhibits some similarities, in particular at high 
temperatures or sufficiently high above the saturation field. The main difference occurs 
at low temperatures and close to the saturation field, where the entropy clearly increases 
from the Heisenberg chain (a) to the J1-J2 chain (b) and further to the sawtooth chain 
(c) - compare also Figs. El and El for more detailed data exactly at He- The increase 
of entropy enhances the magnetocaloric effect as the magnetic field is swept through 
saturation. 

The data in Fig. IHl has been computed numerically for = 20 sites. Hence, it 
is important to discuss finite-size effects. Some related comments for the entropy at 
the saturation field have already been made in the context of Figs. El and We now 
examine finite-size effects for all fields, focusing on the example of the sawtooth chain. 
Fig. [71 corresponds to a part of Fig.lHfc) at low temperatures and fields and shows curves 
of constant entropy for N = 12, 16 and 20. The lowest curves start at S/N = 0.05 and 
increase in steps of 0.05. At low temperatures and close to He = 2 J we can compare with 
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Figure 6. Entropy of A'^ = 20-sitcs lattices as a function of temperature and magnetic 
field for spin-1/2 models: (a) the Heisenberg ciiain, (b) the J1-J2 chain with Ji = J, 
J2 = jJ, and (c) the sawtooth chain with Ji = J, J2 = ^ J. Lines show curves of 
constant entropy starting at the bottom with S/N = 0.05 and increasing in steps of 
0.05. 
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Figure 7. Lines of constant entropy, i.e. adiabatic (de)magnetization curves for tiie 
sawtootli clrain witli Ji = J, J2 = ^ J for N = 12, 16 and 20. 

eq. (j36|) . The corresponding lines of constant entropy are straight hnes since the entropy 
of hard dimers depends on H and T only via the combination {H — H^/T. The lines 
obtained from eq. ()36|) lie on top of the lines obtained numerically for = 20 in Fig. [7| 
in the regime where the latter are linear. Since eq. (j36p is valid for the thermodynamic 
limit, this means that finite-size effects in Fig. [7|are within the width of the lines for 
iV = 20 for H close to Hf.. Further comparison of the curves for different system sizes 
shows that the adiabatic demagnetization curves of a sawtooth chain with = 20 sites 
are essentially free of finite-size effects, except for H < 1.7J, T < 0.15 J. These finite- 
size effects are related to the recovery of a gapless spectrum at T = and for H < 1.5 J. 
Similarly, finite-size effects can also be neglected for the Heisenberg chain and the J1-J2 
chain (Fig. IHfa) and (b)) except for H < 2J, T < 0.15J. Note that in the gapless 
regimes {H < He) temperature varies only weakly with the applied field (compare also 
Fig.lH). 

Finally, Fig. |H1 presents a direct comparison of adiabatic demagnetization curves 
for the spin-1/2 Heisenberg chain, the J1-J2 chain and the sawtooth chain. We start 
with all systems at the same temperature and a magnetic field noticeably above the 
saturation field {H = 3). During adiabatic demagnetization, all three systems cool 
upon approaching the saturation field. Clearly, the J1-J2 chain cools down to lower 
temperatures than the spin-1/2 Heisenberg chain under the same conditions, and the 
sawtooth chain to even lower temperatures. Along with numerical results we show by 
thin lines the analytic results for a free-fermion model for the Heisenberg chain and 
the J1-J2 chain. One observes that free fermions yield an excellent description of the 
Heisenberg chain above the saturation field while slightly bigger deviations for the J1-J2 
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Figure 8. Comparison of adiabatic demagnetization curves when approaching the 
saturation field from above for the spin-1/2 Heisenberg chain, the J1-J2 chain with 
Ji = 1, J2 = 1/4, and for the sawtooth chain with Ji = 1, J2 = 1/2, as indicated in 
the figure. The three bold curves have been computed numerically for TV = 20 sites, 
the two thin curves have been obtained in a free-fermion approximation. 

chain are due to the relevance of interactions in the latter case. 

The adiabatic demagnetization curve in Fig. |H1 of the sawtooth chain corresponds 
to an entropy below the value given in eq. pTj) . Hence, one even observes cooling to 
T— i>0 as if-^ifc = 2. Nevertheless, the example of Fig. |H1 is still at comparably 
high temperatures (entropy). Among others, this ensures that finite-size effects can be 
neglected although the bold curves in Fig. |H1 have been computed with = 20. While 
asymptotic low-temperature expressions like eq. (jlll) are not accurate in this regime, 
the enhancement of the magnetocaloric effect with increasing frustration is evident. 

10. Conclusions and outlook 

When discussing the magnetocaloric effect at low temperatures, one should clearly 
distinguish two different cases: (i) large relative variations of entropy in a magnetic field, 
while the total entropy tends to zero, and (ii) large absolute variations of entropy. Let 
us first discuss the former case, which appears in ordinary nondegenerate quantum spin 
chains. We have shown that an enhanced magnetocaloric effect exists in the vicinity 
of quantum phase transitions in a magnetic field. Specifically, we have studied the 
transition at the saturation field. There are, however, many other systems/transitions 
to which our predictions apply, in particular second-order transitions at the boundaries 
of magnetization plateaux pil21j . Field- induced transitions in ID systems with a singlet 
ground state, like spin ladders or Haldane gap chains, are one specific example (see 
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e.g. |21j for a recent review). Triplet excitations with AS^ = +1 condense, when the 
Zeeman energy exceeds the triplet gap. Triplets can be represented as hard-core bosons, 
which are mapped in one dimension to weakly interacting spinless fermions |12|ll6j. 
Therefore, the entropy of triplet excitations close to the transition field is also described 
by the equations © in the single-particle approximation. A less trivial example is the 
Bose condensation of magnetic excitations in 2D spin systems with ID type degeneracies 
in the energy spectra. Such a possibility occurs in the J1-J2 antiferromagnetic model 
on a square lattice near the saturation field |2ZI as well as in the 2D singlet gap magnet 
Cs3Cr2Br9 |2H|- Triplet excitations above the singlet gap have also a large degeneracy in 
the singlet state of SrCu2 (303)2 j^lHU] which should lead to a large entropy release at 
the gap closing transition. In the frustrated antiferromagnetic chain the magnetocaloric 
effect is enhanced due to a cancellation of the leading term in the one-particle 
dispersion for J2/J1 = 4. Such a cancellation remains valid for ferromagnetic 
sign of the nearest-neighbour exchange Ji < 0. The magnetization curve of the J1-J2 
chain with ferromagnetic Ji < is also particularly steep below the saturation (see e.g. 
Fig. 4 of [201)5 suggesting that the magnetocaloric effect may also be enhanced in this 
regime. Rb2Cu2Mo30i2 is a candidate for a realization in the latter parameter regime 
with a moderate saturation field [3T] . 

A preliminary experimental indication of a large magnetocaloric effect has been 
obtained in pulsed field measurements in a spin-1/2 Heisenberg chain antiferromagnet 
[32] . In the vicinity of the saturation field the magnetization curve deviates significantly 
from the behaviour expected for the corresponding bath temperature indicating a 
sizeable magnetocaloric effect. Thus, the magnetocaloric effect can significantly affect 
experimental results obtained in pulsed magnetic fields. A systematic experimental 
investigation of the magnetocaloric effect near quantum phase transitions can also be 
used with the help of eqs. fj^ . to clarify the role of spin anisotropies in a given 
magnetic material. 

Large absolute variations of entropy in a magnetic field appear for a class of 
geometrically frustrated antiferromagnetic models. The sawtooth chain considered 
in sections [HI and [HI has a finite zero-temperature entropy at the saturation field. 
The analysis of the magnetothermodynamics of this model shows that in a strongly 
correlated spin system there are emergent "paramagnetic" degrees of freedom, which 
experience an effective magnetic field h = {H — H^. Such effective "paramagnetic 
moments" are responsible for an enhanced cooling power in the vicinity of very 
similar to standard paramagnetic salts [2]. However, the "paramagnetic" language is 
not completely adequate and a picture of hard localized magnons is more appropriate for 
the physical description. Among the differences with a simple paramagnetic picture, we 
mention the geometrical origin of the finite entropy: the zero-T entropy at the saturation 
field does not depend on the value of the local spin. Also, there is an asymmetry between 
the /i > and h <Q regions close to the saturation field in Fig. [7| 

The antiferromagnetic sawtooth chain model is realized in delafossite YCUO2.5 jHHj, 
though with a yet unknown ratio of the exchange constants J2I J\- Our theoretical 
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study should also be relevant for 2D and 3D spin models on kagome, pyrochlore 
and garnet (hyper-kagome) lattices, which have numerous experimental realizations. 
The magnetocaloric effect has been known for a long time in gadolinium gallium 
garnet Gd3Ga50i2, which is often utilized in adiabatic demagnetization refrigerators 
|34|l35j. Very recently, a large magnetocaloric effect has also been measured in the 
pyrochlore antiferromagnet Gd2Ti207 jHE]- Magnetic Gd^"*" ions have a rather large 
spin S = 7/2. This justifies the classical approximation applied for pyrochlore and 
garnet antiferromagnets in jH]. Still, the classical model cannot describe the behaviour 
at low temperatures T < JS. One of the physical implications of the present work 
is that the entropy on frustrated lattices has a purely geometrical origin. Therefore, 
5* = 1/2 frustrated magnetic materials should be as good as S ^ 1/2 compounds for 
cooling by adiabatic demagnetization, contrary to the experience with paramagnetic 
salts j2]. We suggest that a synthesis of new spin-1/2 antiferromagnetic materials with 
kagome, pyrochlore and garnet lattices should be interesting both from a fundamental 
point of view and with regard to possible technological applications in the field of low- 
temperature magnetic cooling. 
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